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Abstract 

We investigate shock formation in an asymmetric exclusion process with creation and annihi- 
lation of particles in the bulk. We show how the continuum mean-field equations can be studied 
analytically and hence derive the phase diagrams of the model. In the large system-size limit direct 
simulations of the model show that the stationary state is correctly described by the mean-field 
equations, thus the predicted mean field phase diagrams are expected to be exact. The emergence 
of shocks and the structure of the phase diagram are discussed. We also analyse the fluctuations 
of the shock position by using a phenomenological random walk picture of the shock dynamics. 
The stationary distribution of shock positions is calculated, by virtue of which the numerically 
determined finite-size scaling behaviour of the shock width is explained. 



I. INTRODUCTION 



The analysis of self-driven particle-models is a central issue of non-equilibrium statistical 
mechanics 11, These models show a variety of generic non-equilibrium effects, 

in particular in low dimensions. Although the behaviour is generally rather complex, some 
models are simple enough to be analysed in great detail. In rare cases it is even possible 
to obtain exact results for the stationary state . A very prominent example of such a 
model is the asymmetric exclusion process (ASEP) which comprises particles hopping in a 
preferred direction under the constraint that they cannot occupy the same lattice site. Exac ; 
solutions for the stationary state exist for periodic lUfJ and open boundary conditions 
and different update schemes 

BHQ. TWefo^odel^ beeped in o ri eT t o 
develop more general concepts for systems far from equilibrium, e.g. a free-energy formalism 
[l^ . Variants of the model have also allowed the study of shocks which are discontinuities 
in the density of particles over a microscopic distance 13| . 

The ASEP and related models are not only of academic interest, but have a number 
of important applications, e.g. as simplified traffic models Q. Here we are interested in 
a variant of the ASEP, which is motivated by biological transport processes in living cell 
systems, where particle non-conservation in the bulk of the system is allowed [3| (see also 

An important feature of living cells is their ability to move and to generate forces Il5l 
Hi- On a microscopic level these forces are mostly generated by motor-proteins 

HQ, 

which are able to perform directed motion along one-dimensional paths or filaments. 

Many different motor proteins can be distinguished [rsj |. but they have a few common 
properties: a head, which can couple to a filament, where it performs a directed stochastic 
motion; and a tail, which is attached to a specific load, which has to be transported through 
the cell. The coupling of the motor protein heads to the filament is reversible, thus the 
motor proteins will attach to the filament, perform stochastic directed motion for some time 
and eventually detach from the filament. The typical distance between attachment and 
detachment of the motor protein to the filament depends on the particular type of filament. 

To a first approximation the motion of many motor proteins along a filament can be 
modelled by the asymmetric exclusion process. An important feature, which is not described 
by the ASEP, is the attachment and detachment of the motor-protein heads. This feature 
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has been included in a recent model by Parmeggiani et al (see also |28j), which can 
be viewed as a grand-canonical counterpart of the ASEP in the sense that in the bulk the 
particle number is not conserved. 

For open boundary conditions, however, the situation is different. In this case, if the rates 
of attachment scale correctly with system length, one observes a subtle interplay between 
the left, right and bulk particle reservoirs. If we fix the attachment and detachment rates 
then there are whole regions of the phase diagram (spanned by the densities of the boundary 
reservoirs) where one observes the localisation of shocks in the bulk of the system. This is in 
contrast to the ASEP, where shocks move with constant velocity and are generally driven to 
the boundary of an open system. The shock has zero velocity only on the phase boundary 
where two phases of different density co-exist. Even in this case the shock is not localised 
since it diffusively explores the whole system. 

The physical origin of the shock localisation in the present model and a discussion of 
the phase diagram, in particular the phase where shocks appear, are the subjects of this 
article. We shall show how the phase diagram can be predicted through simple considerations 
pertaining to a continuum mean field description which only retains first order terms i.e. the 
phase diagram can be predicted through the study of a simple first-order nonlinear partial 
differential equation. Going beyond the mean level, in the shock phase we describe the 
dynamics of the shock by a random walker with space-dependent hopping rates. In this way 
the localisation of the shock can be understood. 

The article is organised as follows. In the next section we give the definition of the model 
and introduce the stationary solution in case of periodic boundary conditions. In section ITTT1 
we introduce the mean field equations for the open system and discuss their solution by 
means of characteristics. Then we discuss the stationary solutions on the mean field level 
and compare the mean field results to Monte Carlo simulations. Fluctuations of the shock 
positions are analysed in section and finally some concluding remarks are given in the last 
section. 

II. MODEL 

We consider a one-dimensional open chain of N sites, which can either be empty (ji = 0) 
or be occupied with one particle (t, = 1). Particles can jump to the neighbouring site if 
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it is empty. In addition, the bulk sites are coupled to particle reservoirs, i.e. particles are 
attached with rate uja and deleted with rate ujd- The particle reservoirs at the boundaries 
are different from the bulk. At the first site particles are attached with rate a and deleted 
with rate Ur>, while at the last site particles attach with rate u A and are deleted with rate 
(3 jl^j]. By a rate r it is meant that in infinitesimal time interval dt the probability of the 
event occurring is rdt. Schematically the dynamics can be written as follows: 

*0 ^ *1 (1) 
1*^0* (2) 
10 -^U 01 (3) 

where 1(0) corresponds to an empty (occupied) site and * implies that the update is inde- 
pendent of the state. 

The dynamics at the left hand boundary (site 1) is insertion of particles 

0- ^1 (4) 
and at the right hand boundary (site N) is removal of particles 

1- ^0. (5) 

The exact equation for the evolution of the particle densities (rj) away from the boundaries 
(1 < i < N) is given by 

^ = (r f _i(i - Ti)) - (n(l - 7-i+O) + u A (i - n ) - uo D (n) , (6) 

where (...) denotes the statistical average. At the boundaries the densities evolve as: 



dt 

d(r N ) 
dt 



-(n(l - r 2 )> + a(l - n) - LU D (n) , (7) 
(TN-iil - t n )) + u A (l - t n ) -P{t n ) . (8) 



First consider the steady state on a periodic system (where site A^+l is identified with site 
1 and (|4|5jl do not apply). Assuming translational invariance, © is satisfied by (rj) = p cq 
where 

Peq = U A j (U A + UJ D ) . (9) 
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We refer to the density as the equilibrium density as it is the density obtained in the 
Langmuir absorption model j^. Furthermore, it can be verified that the steady state of 
the system is given by a product measure with density of particles p eq . 

We can apply particle-hole symmetry in order to simplify the discussion of the model. In 
the case of open boundary conditions, the system of equations (JOHHJ) is invariant under the 
simultaneous exchanges a ^ j3, u>a ^ i ^ L — i and pi «-> 1 — pi. 



III. MEAN FIELD EQUATIONS AND CHARACTERISTICS 

In the large N limit we can make the continuum mean field approximation to (Kfl). First 
we factorise correlation functions by replacing (rj(l — r i+1 )) with (rf_i)(l — (Tj}) [21[ then 
we set 

1 dp 1 d 2 p , , 

<^> = ^ ± j?S + 2jPa?-- (10) 

Keeping leading order terms in 1/N, one obtains 

^ = -(1 - 2p)|| + u D N [K - (1 + K)p\ , (11) 

where r = t/N and 

if = uj A /uj D . (12) 
For the open system we shall be interested in the scaling limit where 

Qa — ^aN and Qe> = uj£>N (13) 

are finite as iV —>■ oo. The boundary conditions become p(x = 0) = a and p(x = 1) = (1—0). 
To understand the first-order differential equation (jllj) one can study the characteristics 



221 ] . which are defined for a quasi-linear equation, 



r,p)^ + b(x, r, p)|^ = c(x, r, p) , (14) 



by the equations 

dx b(x, t, p) dp c(x, r, p) 



(15) 



dr a(x, t, p) dr a(x, r, p) 
Roughly speaking, characteristics are curves along which information about the solution 
propagates from the boundary conditions of the partial differential equation. Based on 



characteristics, Lighthill and Whitham |23J] developed the theory of kinematic waves for 
equations with mass conservation and showed how kinematic shock waves arise. Here we 
generalise this picture to equation where the number of particles is not conserved. 
In the present case the characteristics are given by 

1 = 1-2, (16) 
|£ = Si D (AT - (1 + K)p) . (17) 

I 1 

Equations (j 1 7j) are to be interpreted as kinematic waves [23], which propagate changes 
in the density, moving with speed 1 — 2p but with the density of the wave p itself changing 
with time. In the absence of creation and annihilation of particles the density of the wave 
is constant in time and the waves propagate in straight lines. However in the presence of 
creation and annihilation of particles the waves will follow a curve in the x—t plane. For 
example, consider a density fluctuation starting at the left boundary with p(0) = a < 1/2 
and a < K/(l + K). Initially the fluctuation will propagate to the right with speed 1 — 2p 
with its density increasing and speed decreasing. If K/(l+K) < 1/2 i.e. K < 1 the density 
will approach p = K/(l + K) and the fluctuation will propagate to the right with a fixed 
speed. However if K/(l + K) > 1/2, after some time the density of the fluctuation will reach 
p = 1/2 and the fluctuation will cease to propagate. Similarly, a kinematic wave starting at 
the right boundary with p(l) = 1-/3 where (3 < 1/2 and (3 < 1/(1 + K) will travel to the 
left with decreasing density and decreasing speed. 

When two characteristic lines cross, multivalued densities are implied, therefore the de- 
scription by a first order differential equation breaks down. However Lighthill and 
Whitham showed that the effect is that a shock, i.e. a discontinuity between the densities pi 
and p2, arises at the meeting points of the two characteristics and this discontinuity travels 
with speed v s . This speed is determined by balance of mass current to be [23 

p 2 (l-p 2 )-Pi(l-Pi) 1 , 1H x 
v s = = 1 - Pi - Pi • (18) 

P2 ~ Pi 

In the present case, although the mass is not conserved, the mass current between through 
the shock still implies that its velocity is given by (|18j) . Thus for a<l/2, /5<l/2we have 
the possibility of a shock forming then being driven to a position where the mass current 
through it is zero and the shock remains stationary. In the next section we shall show how 
this picture is borne out by solving the steady state mean-field equation. 



IV. STEADY STATE SOLUTION 



Setting the time derivative of (fTTj) to zero yields 

dp 



[l-2p)^--n D [K-(1 + K)p]=0 



(19) 



This is a first order ordinary differential equation, which in principle can be solved ana- 
lytically. The only difficulty is the occurrence of shocks in the solution. To construct the 
solution we integrate from the left boundary (p(0) = a) to find a profile pi{x): 

1 - 2p 



x 



I rpi( x ) 

?W Q dP K-(l + K)p 

K 



n D (i + K) 



2{pi -a) + 



In 



K 



K) P i 



K - (1 + K)a 



and integrate from the right boundary (p(l) = 1 — (3) to find a profile p T (x) through 



1 



x 



K - 1 

2(l-/3-p r ) + - — In 



K-(l+K)(l-p) 



K-(1 + K)p r 



(20) 
(21) 

(22) 



To determine the full profile across the system we have generally to match these two profiles 
at a shock whose position is to be determined. 



The case K = 1 



First we consider the special case of K — 1, i.e. where the attachment and detachment 
processes have equal rates. The steady-state mean-field equation (fEU) reads 



(2 P -i)(g-n) = o, 



(23) 



where Q = Qa = This can be solved explicitly by a piecewise linear trial function 

leading to the condition: either p(x) = const = p eq = | or p(x) has slope Q. Note that this 
solution implies that higher derivatives in eq. (jl())l vanish exactly. 

First we consider the parameter regime a < 1/2, (3 < 1/2, where shock formation is 
possible. In order for the shock's position to be stable i.e. for it not to be driven out of the 
system the shock's speed, as determined by (fTHj) . should be zero. Thus the densities at the 
discontinuity should be related by p r (x s ) = 1 — pi(x s ) and this determines the position of 
the shock x s : 



X s 



13 -a 1 
2VL + 2' 



(24) 
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FIG. 1: Density (r<) (left) and flow (#) = {n(l - r i+1 )) (right) profiles for a = 0.2, /? = 0.3, 
K = 1, and O = 0.1. For the MC results we set x = i/N . The continuum mean-field theory results 
are compared to Monte Carlo (MC) simulations of different system sizes. The agreement between 
MF and MC results is improved for larger system sizes. 
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FIG. 2: Same as fig. H but the parameters a = 0.2, (3 = 0.1, K = 1 and O = 2.0 are used. For this 
larger value of O = 2.0 shock formation is preempted by reaching the density p eq = 1/2 see ©• 



The height of the shock A is given by 

A = Pr (x s ) - P i(x s ) = i-(a + p)-n = Q c -n 



(25) 



The last equation can be used in order to discuss the parameter dependence of the model. 
If A is positive and < x s < 1, we find indeed a shock, which connects two domains with 
linear density dependence, as illustrated in fig. ^ For Q > Q c one does not observe shock 



formation, but a section 



l-2a 



< X < 1 



1-2/3 



where p 



i 



,scc 



fig. m ■ 



In the ASEP the line a = ft is a phase boundary where shocks, between a high density 



S 



region coexisting with a low density region, exist: Mean field theory predicts the shock is at 
x = 1/2 although the exact solution shows that the shock is actually delocalised and yields 
a linear density profile Q- In the present case, although the shock's position is x = 1/2 
when a = (3, the role of this line changes and one does not observe a phase transition in 
crossing it. 

Moreover linear density profiles, which are observed for a + (3 + Q = 1, do not signify 
phase coexistence as for the ASEP, but indicate a vanishing height of the shock (A = 0). 
Also note that a strictly linear profile is observed only in the limit N — > oo in contrast to 
the the disorder line of the ASEP (which is recovered in case of Q = 0) where the constant 
density solution is valid for finite systems as well. 
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FIG. 3: Same as fig. ^ but the parameters a = 0.8, (3 = 0.2, K = 1 and SI = 0.6 are used. The 
average densities at the first site are indicated by symbols. 



Next we discuss the case a > 1/2 and (3 < 1/2. In this parameter regime the density 
profile is given by 



p[x) 



max 



for x > 0. 



(26) 



Equation (j2T?j) implies a singularity at x = since the boundary condition is p(0) = a. For 
finite systems the singularity at x = is softened as shown in figure fig. El The behaviour 
for a > 1/2 and (3 < 1/2 can be obtained from particle-hole symmetry. 
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Finally in the maximal current regime the density profile in the limit L — > oo is given by 



a for x = 
PO) = I 1/2 for < ac < 1 
1 — (3 for x = 1. 

i.e., apart from the two singular points x = and x — 1, the density is constant. 



(27) 




FIG. 4: Phase-diagram for K = 1 and < 1/2. Indicated are the high- (H) and low-density (L) 
phase, the shock phase (S), the maximal current phase (M), and coexistence between a maximal 
current and low- (LM) or high-density (HM) domain. The LMH-phase indicates the presence of 
three domains. A linear density profile is observed at the transition line between the LMH and S 
phase. 



The above results lead to the phase diagram shown in fig. Three phases can be 
distinguished, where the density profile does not reach the equilibrium p eq : The high and 
low density phase, where only a single domain exists in the system, the shock (see fig. Q 
for the corresponding density profile) region, where high and low density domain coexist. 
Then there are four phases, where the density profile has a section of constant density: 
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The maximal current phase (M) is obtained for a > 1/2 and (3 > 1/2. If a > 1/2 and 
1/2 > /? > 1/2 — Q we observe coexistence between a section of density 1/2 at the left and 
a linear profile as shown in fig. El This phase is indicated by (HM) in the phase diagram, 
the corresponding low density phase by (LM). Finally a section of constant density may 
coexist with a high- and low-density section (LMH) (see fig. |2J) • We observe all phases only 
if Q < 1/2. For l/2<fi<lthe high- and low density phases vanish, while for Q > 1 the 
formation of a shock is excluded. Our MC analysis illustrates the validity of the mean field 
results for large system. This is in contrast to the case Q = 0, where e.g. for a = j3 the 
exact density profile even in the limit iV — > oo is different from the mean field results. 

B. The case K / 1 

In the case K ^ 1 the equilibrium density is different from p = 1/2. Therefore it is 
impossible to observe a maximal current phase, because the bulk particle reservoirs destroy 
any maximal-current domain. The absence of the maximal current phase originates in 
(I16I17J) : For the kinematic wave to be stationary one requires both (|16I17|) to be zero i.e. 
a bulk density satisfying both p = 1/2 and p = It is also important to notice that 

the solution of the mean field equations is not piecewise linear, i.e. higher order terms of 
eq. (jl()j) do not vanish. We now analyse the mean field equations for K > 1 in more detail. 
The corresponding results for K < 1 can be obtained from particle-hole symmetry. We 
consider separately the regime where a < 1/2 and /3 < 1/2 and the complementary regime. 

Case a < 1/2 and (3 < 1/2 

For low values of a and f3 we expect the existence of a high-and low-density phase as well 
as the formation of shocks within a certain density regime as for the special case K — 1. 
The transition lines can be obtained by analysing the shock position, which is determined 
from the condition p r = 1 — p\. The shock separates a region where the density is given by 
(121 J) and a region where the density is given by ()22j). 

For low a the solution ()21|) may propagate all the way to the right boundary. If the 
density at the right boundary satisfies pi(l) < Pr(l) = P any shock will be driven out and 
the system will be dominated by pi. This will be referred to as a low density region (LD). 
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FIG. 5: Density profiles obtained from Monte Carlo simulations of the model for a = 0.2, (3 = 0.6, 
K = 3, £Id = 0.1 and N = 100, 1000, 10000. The full line gives density profiles obtained from the 
continuum mean field theory, equations (|21|) and (|22|) with (3 set equal to 1/2. 



The transition line between LD and S regions is given by the condition pi(l) = f3: 



2{(3 - a) + 



K - 1 



In 



K-(l + K)/3 



K-(1 + K)a 



(28) 



Similarly for large (3 the solution p r may propagate all the way to the left boundary where 
its density is p r (0). A shock between p r and pi will be driven to the left hand boundary if 
p r (0) > 1 — pi(0) = 1 — a. This will be referred to as the high density region (HD). The 
transition line between HD and S regions is given by the condition p r (0) = 1 — a: 



n D (l + K) 



K-(l + K)(l-/3) 



K 



;i + #)(!-«) 



(29) 



Cases a > 1/2 or (3 > 1/2 

For a > 1/2 or (3 > 1/2 second-order terms (e.g. terms involving d 2 p/dx 2 ) have to be 
retained in (jll)) for steady state solutions to exist. However the effect is, for example in the 
case a < 1/2,(3 > 1/2 at the right hand boundary, that over a finite distance (order 1/N in 
terms of the variable x) the second order terms will match the density with that implied by 
p r . Thus the effect is that there is no shock and (3 > 1/2 may for the purposes of the phase 
diagram may be effectively considered as (3 = 1/2. In fig. El comparison of direct simulations 
of the model when (3 > 1/2 with density profiles obtained from the continuum theory (with 
(3 effectively considered as 1/2) indicate the validity of our approach for in the large iV limit. 
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FIG. 6: Mean- field phase diagram for O/j = 0.1 K = 3. (right). The phase boundaries between 
the low density (L), shock (S) and high density (H) phases are calculated by using l(25|) and 
Our results are in excellent agreement with the findings of ^| who analysed numerically the 
second-order mean field equations. 



For finite N, finite size effects can partially be included by considering second-order 
terms in the mean-field description. Deviations between second-order mean-field and the 
exact results are due to the fluctuations of the shock, which are not correctly described by 
mean-field theory and have to be treated separately. This will be done in the following 
section. However at this stage we note that the shock is indeed localised and that the width 
of the shock growth subextensively, i.e. the shock is sharp in the limit iV — > oo. 



Phase diagram 

Our analysis for K ^ 1 leads to the phase diagram shown in fig. El We distinguish only 
three phases, i.e. the high- and low-density phase as well as a the formation of shocks in 
an intermediate parameter regime. Compared to the case K — 1 the number of phases is 
considerably reduced, due to the absence of a maximal current phase as discussed above. 

As the shock is sharp the fluctuations of the shock, which we treat in the next section, 
do not affect the phase boundaries. Therefore we expect that the true phase diagram of the 
model is represented by fig. H3 
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V. DYNAMICS OF THE SHOCK 



A qualitative understanding of the shock dynamics can be easily obtained from mass 
conservation, i.e. by means of the continuity equation. Compared to the ASEP additional 
source and sink terms have to be introduced, reflecting the on-site input and output of 



24j . who generalised the 



particles. This has been done in a recent work by Popkov et al. 
domain wall picture to models with particle in- and output, where the continuity equation 
is given by: 

f)n( nr f 1 f) 

n A (i-p(x,t))-n DP (x,t) . (30) 



dp(x,t) d 



By inserting the known relation j(p) = p (1 — p) they recovered eq. (jlljl . The importance of 
the analysis by means of the hydrodynamic equation is given by the fact, that this approach 
can be used for general models if the flow-density relation is known 24 1. 
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FIG. 7: Left: Density profiles in comparison with the DW predictions. The chosen parameters are 
a = 0.2, (3 = 0.3, £1 = 0.2 and K = 1. Right: Flow profile of the left (solid line) and right domain 
(dashed line) for same set of parameters. The arrows indicate the bias of the random walker. 

A step beyond the mean field level is provided by interpreting the shock as a random 
walker. For the ASEP on the line a = (3 the shock dynamics can be modelled as an unbiased 
random walk j^j and this allows one to calculate its diffusion constant. Moreover one can 

26, Ez|. In 



generally consider kinematic shock wave dynamics as a biased random walker 
the present case we model the shock by a random walker with the site-dependent hopping 
rates 

j-( x ) / • \ J+( x ) 



Wi[Z 



AW 



Writ) 



AW 



(31) 
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where wi (w r ) denotes the hopping rate to the left (right), the flux in the low 

(high) density domain at site i and A(i) the height of the shock at position i. Modelling the 
shock dynamics as a random walk with reflecting boundary conditions allows us to derive 
analytical expressions for the stationary distribution p s (i) of shock positions. The stationary 
distribution p s (i) has to fulfil the condition 

w r {i)p s {i) = wi(i + l)p s (i + 1). (32) 

One can solve this discrete equation explicitly, but it is simplest to proceed by making a 
continuum approximation. We expand the probability distribution to first order in 1/N, as 
in eq. (|1U|). and use the stationarity condition to obtain the differential equation: 

y'{x) + Ny(x) (l - ^l) = , (33) 



(34) 



Wi{X) 

where y(x) = p(x)wi(x). The solution of this equation is given by 



p(x) 1 



-N I [ 1 - 1 dx' 



x V wi^r 



where M = J Q p(x)dx. Explicit expressions of the distribution can be given in case of K = 1 
which will be discussed in detail. For this case we get the unnormalised distribution 

iV(l-A)A 1 v JV ' 1 + A > A -1 

« ,) = ( ,+ 5) ' (t"') ~e~^\ (35) 

where C = 4NQ ( 1 _ A ^ 1+A - ) ■ The Gaussian is obtained from a logarithmic expansion of the 
exponent, which is justified in the limit N — > oo j^]. 

Using the random walk picture, the density profile can be obtained from the shock dis- 
tribution p(x) through 

/•X pi 

p{x) — p r (x) / p{x')dx' + pi(x) / p(x')dx'. (36) 

JO Jx 

The Gaussian approximation of p(x) leads to the following form of the density profile 

+ tlx + a . (37) 




NQA 

erf|2W ;i + A)(i-A) (x -^ 



In order to check the validity of this simple phenomenological picture we compare the an- 
alytic predictions (|37|) for different density profiles in the (S) phase with results of MC 
simulations [30] • Fig. shows that the DW predictions are in good agreement with the MC 
results, and the accuracy of the DW description improves for larger system sizes. 



From (JBTj) it jo 
numerically in 



lows that the width of the shock scales as iV v with v — 1/2 as was found 
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VI. CONCLUSION 

The model studied by Parmeggiani, Franosch and Frey |14| can be interpreted as a generic 
model for the collective behaviour of molecular motors. The physics of the model is governed 
by the competition of different particle reservoirs and the stationary flow of the self-driven 
particles. In the case of periodic boundary conditions, one easily verifies that the stationary 
state of the process is described by a product measure. More interesting features are observed 
in the case of open boundary conditions. Here, when the rates for attachment and deletion of 
particles in the bulk are appropriately scaled with system size, localisation of a shock arises 
between the region of the system controlled by the left boundary and the region controlled 
by the right boundary. 

Most of the features of this model can be explained by a mean-field analysis, which we 
believe to be correct in the limit of large system size. This is supported by figs 1,2,3,5 where 
direct simulations for the density profiles converge for large system sizes to our mean-field 
predictions. In view of this we believe the mean-field phase diagrams figs 4,6 are in fact 
exact. 

By considering the characteristics of the mean field equations the formation and locali- 
sation of the shock can be explained: The characteristic solutions propagating from the left 
and right boundaries are matched at the shock whose position is fixed by the condition that 
the mass current through the shock is zero. In the presence of a shock, the leading finite- 
size corrections are due to the fluctuations of the shock position, which can be described by 
mapping the dynamics of the shock to a random walk with site dependent hopping rates. 

Apart from the importance of the ASEP with creation and annihilation of particles in 
the bulk as a generic model for molecular motors it is of special interest for the general 
formalism of non-equilibrium statistical mechanics. In particular one can inter pret the bulk 
non- conservation of particles as exchange of particles with a bulk reservoir 28|. In this 
way the model can be thought of as a grand-canonical counterpart to the ASEP. As the 
ASEP and its variants can be analysed in detail, they might help us to understand aspects 
of different ensembles in the context of non-equilibrium statistical mechanics. 
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